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Abstract 

Wavelets are an exciting new topic in applied mathematics and 
signal processing. This paper will provide a brief review of wavelets 
which are also known as families of functions with an emphasis on 
interpretation rather than rigor. We will devise an indirect use of 
wavelets for the solution of integral equations based upon techniques 
adapted from image processing. Examples for resistive strips will be 
given illustrating the effect of these techniques as well as there promise 
in reducing dramatically the storage requirement in order to solve an 
integral equation for large bodies. We also will present a direct im- 
plementation of wavelets to solve an integral equation. Both methods 
suggest future research topics and may hold promise for a variety of 
uses in computational electromagnetics. 

1 Introduction 

Wavelets have generated significant excitement among applied mathemati- 
cians and engineers recently due to their unique properties and diversity. 
Wavelets are actually a family of basis functions which are generated by 
translation and dilation of a single function. The uses are diverse since many 
different wavelets may be generated which possess some properties which 
are particularly well suited for a given purpose. Three prominent uses of 
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wavelets are image compression which is discussed by Mallat( [ 1] ,[2] ) and 
Daubechies[3], solution of partial differential equations by the staff at Aware 
Inc. ([4], [5]), and the solution of Fredholm Integral Equations by several 
members of the Yale research staff ([6], [7], [8]). 

This paper will give a brief and basic review of wavelets with an empha- 
sis on the character of these functions rather than a rigorous mathematical 
treatment. Our purpose is to apply wavelet expansions or wavelet trans- 
forms to the solution of electromagnetic problems. We will show that a novel 
technique which is adopted from image processing may be applied to tra- 
ditional Method of Moments(MoM) computer codes to achieve significant 
storage savings. Since modern computers are becoming increasingly fast, the 
primary limitation on MoM implementations is the 0{N 2 ) storage require- 
ment. The use of wavelets is intended to alleviate this requirement. This 
“image processing’’ technique may be used with existing MoM implementa- 
tions to reduce storage. Finally we will look at direct applications of wavelets 
in electromagnetics with a modal solution of the Fredholm equation and dis- 
cuss future research topics. Although we have only used E-polarized resistive 
strip examples for this paper these techniques are applicable to other elec- 
tromagnetic problems of interest. In particular, although we only use MoM 
examples herein, we are interested in applying wavelet analysis in the solution 
of the boundary integral associated with the FEM-BE method. 


2 Wavelets and Dilation Equations 

Wavelets and their underlying form, a dilation equation, is an unfamiliar 
concept to most engineers. We will attempt to give a feel for what a wavelet 
is and what properties it possesses which are attractive for electromagnetics 
without going into detail. Strang[9] gave an excellent high-level review of 
the properties of wavelets while DaubechiesflO] provides a rigorous review 
with significant new material. Mallat[2] was able to relate wavelet decompo- 
sition and reconstruction in terms of multiresolution analysis which allowed 
significant progress to be made in exploiting wavelets for the aforementioned 
applications. 

In order to understand wavelets, we must first look at a dilation equation 

2x-k) ( 1 ) 
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which is recognized as a two-scale difference equation. Strang[9] reviews 
the criterion which must be met in order for (1) to be unique. The scal- 
ing function”, is determined recursively and thus many properties of 
it must be observed rather than defined. Daubechiesjl 0] was able to con- 
struct orthonormal wavelets with compact support m both domains which 
was thought to be impossible by using the recursion 

<j)j{ x) = y; c k <f>j- i(2j — k) (2) 

k 

with the box function as the fundamental function ( 0 D { ? )) and certain simple 
conditions which determine c^. The family of w’avelets derived b\ Daubechies 
will be used for the examples of this paper although this is certainly not the 
only or necessarily optimal choice of w’avelets. 

At this point, w T e introduce the continuous wavelet family w’hich comes 

from the scaling function, 4>{x), as 

ll'(x) = y/s<f>(s{x - ti) ) (3) 

where the scale parameter(s) and the translation parameter(u) range over the 
positive real axis and the entire real axis, respectively (The notation adopted 
for the remainder of this paper is from Mallat [1]). A fundamental property 
of wavelets is their approximation capability. Daubechies family is denoted 
by the number of recursions used to determine the analyzing wavelet. For 
example we will denote the Daubechies w'avelet which has four coefficients 
as DAUB4 which can recover up to a linear polynomial since the two low r est- 
order moment vanishes as discussed by Strang[9]. As terms are increased, 
more moments vanish and thus the polynomial approximation increases in 
order (e.g. DAUB6 will recover a quadratic). 

As was previously mentioned, the Daubechies family is not only orthonor- 
mal but they enjoy the useful property of compact support. For example, 
figure 1 illustrates DAUB20 wavelet compared with the comparable cosine 
if we wish to represent an impulse “frequency” of / = y^ 4 a w here A is the 
spatial sample size. As seen, the cosine( Fourier series basis) is compact in 
the “frequency” domain yet requires infinite support in the spatial domain. 
However, the Daubechies w’avelet is compact in both domains. A further 
example is given in figure 2 where the excitation impulse is at / = 
and we observe that as the “frequency” increases, the analyzing wavelet has 
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diminishing spatial support. This is the most important property of wavelets 
since finer spatial resolution is used for “high frequency” components which 
results in a more accurate evaluation than is possible for the comparable 
Fourier coefficient which is computed with constant sample size regardless of 
the “frequency”. 

Several groups ([4], [5], [6], [7], [8]) have recognized that the property of 
compact support in both domains will result in a sparse matrix when an 
integral operator is converted to a matrix operator. We will discuss in the 
next section the fact that if one thinks in terms of multiresolution analysis, 
the preceding observation is obvious. 

3 Matrix Compression Techniques 

One of the most promising applications of wavelets is image compression 
for the transmission and storage of detailed images. For example, The FBI 
is looking at using wavelet compression techniques in order to store their 
extensive finger print library without loss of fidelity. We will give a brief 
review of the wavelet transform, multiresolution analysis , and what this 
topic has to do with electromagnetics. 

The continuous wavelet transform is given here without proof[l] 

F(s,u) = f f(x)r(>,(x - u)dx (4) 

where the wavelet family is given 

tMx) = Vsrp(sx) (5) 

and i>(x) is the mother wavelet which is denoted as <f>(x) in (8). In order to 
reconstruct a signal from its wavelet transform, we define the inverse wavelet 
transform as 

/ oo roc 

/ F(s,u)xl>.(x - u)dsdu (6) 

-oo JO 

A more practical transform for numerical applications is the Discrete Wavelet 
Transform (DWT) which acts on a sequence of samples. In order to cover the 
phase-space, we need to uniformly sample the translation parameter while 
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exponentially sampling the scale parameter[l]. The resulting DWT is given 


roo 

Wd[/(j>)] = / J(t)v v (u - na J )dx 

J -oo 


(<) 


where a is the elementary dilation step. The signal may he reconstructed 
from its DWT by the series 

f(x) = f( u )^'v( u - riQ_J ) > t' 2 j(u - rio' J ) (8) 

J n 

with < •,• > denoting the standard inner product. In practice, we choose 
a = 2 for historical reasons although this is not necessary. 

The Fourier transform is recognized as a transformation from the spatial 
domain to the wavenumber domain. However, the wavelet transform is a 
transformation from the spatial domain to the scale domain. This is the heart 
of the multiresolution interpretation of (7) as well as its fast implementation, 
the Fast Wavelet Transform(FWT). Each successive coefficient of the DWT 
describes the difference or details between the current level and a the previous 
decimated level. Figure 3 illustrates this concept. The first level, figure 
3a, describes the smooth component of the signal. Each successive level 
adds detail to the reconstruction. Thus a smooth signal will have a large 
large-scale component and very small detail components which results in 
signal compaction. The efficiency of the FWT algorithm follows from the 
pyramid structure shown in figure 3. The highest detailed level (smallest- 
scale) requires that all samples be used in its computation. However each 
successive level requires half (a = 2) the number of samples due to dilation 
which results in the FWT having the same order of operations as the FFT. 

Mallat[2] was the first to recognize that this property naturally satis- 
fies the procedure used in pyramid-encoding schemes. If one applies a two- 
dimensional FWT which is computed ip a manner similar to the most com- 
mon implementation of the two-dimensional FFT to an image, the image 
can be compressed into few large components and many components which 
are less than some tolerance whose omission does not adversely affect the 
fidelity of the image. Significant compression can occur depending on which 
analyzing wavelet is used. 

We realized that the typical impedance matrix generated by a MoM code 
may be thought of as a complex-valued image! Thus it is reasonable to as- 
sume that we may apply the wavelet transform to a dense impedance matrix 
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to generate a sparse matrix. It should be mentioned that this technique 
is similar in concept to the Impedance Matrix Localization! IML) method 
proposed by Canningfll] who uses physical arguments to generate a trans- 
formation matrix which yields a sparse matrix by using sub-domain basis 
functions. It is not yet apparent which method yield superior results. 

Suppose one wants to solve for the E-polanzed scattering from a resistive 
strip. The appropriate integral equation is a Fredholm equation of the Second 
Kind 

= R{x,y)J 2 {x,y) + J J c M* *y)Hi l) [k 0 \y/[x - x') 2 + (y - y'|) 2 ]d/' (9) 

This may be solved numerically by converting the integral operator into a 
matrix operator via MoM to form the matrix equation 

Z -A = E l (10) 

where each entry of J z and E' z corresponds to a spatial sample. If we look at 
the impedance matrix, Z, typically it is dense. This is shown in figure 4 which 
is the impedance matrix generated by pulse basis-point matching on a lAxlA 
corner reflector where clearly all elements of this matrix are important. 

If however we apply a wavelet transform to both sides of (10), we get the 
transformed equation in the scale-domain 


ZJ. = Z (11) 

Figure 5 illustrates the matrix shown in figure 4 after wavelet transformation 
by the DAUB10 analyzing wavelet. Clearly, this matrix has fewer signifi- 
cant entries and is thus a sparse matrix. For images, edges tend to be very 
important thus we would use very compact (low-order) wavelets in order pre- 
serve these edges. However, when one has a smooth matrix such as the one 
shown in figure 4, smooth wavelets (large-order) will tend to have greater 
compaction capability than a small-order wavelet. In order to verify this 
claim, we looked at the compression realized by using different wavelets from 
the Daubechies family and a truncation tolerance of 10 -2 . Table 1 summa- 
rizes the results of this study for the corner reflector when the matrix size 
was 32x32. We observe that the apparent storage was reduced to C?(20N) 
from jV 2 without loss of fidelity as shown in figure 6. However, the large- 
order wavelets did not exhibit a significant advantage with respect to the 
small-order wavelets. 


6 



The reason for this unexpected result is the fact that the matrix in ques- 
tion is not large. A similar comparison is shown in table 2 for the case 
of a 6A curved(inclusion angle = 45°) metallic strip which used a 128x128 
impedance matrix. We recognize some advantage using large-order wavelets 
while maintaining fidelity as seen in figure 7. As a final example, we look 
at the simulation of a metallic half plane-linearly tapered resistive half plane 
junction which used a 1024x1024 matrix. A clear advantage using large-order 
wavelets is seen in table 3 where the apparent storage is reduced to 0(7ON) 
which results in a 93% storage savings without loss of fidelity as shown in 

figure 8! .... 

The astute engineer will recognize the use of the term “apparent in the 

preceding storage analysis. This is due to the fact that you need to generate 
the entire impedance matrix before transformation which results in G{ A ) 
storage. However, we have devised a simple block-by-block technique which 
will result in an actual memory savings. If one writes (10) as an augment ed 
matrix equation where there are (TV/ A') 2 blocks of A' 2 elements each where 
K is of course a power of two(a = 2) then we may write (11) as 

N/K 

= ft / = l-.A'/A' (12) 

k=l 

To illustrate this technique, we looked at the junction simulation previously 
given. Table 4 compares the compression results realized when K - 64 and 
the tolerance was set at 10~ 4 . The total storage requirement is A' 2 + 0(240 A' ) 
which is a 75% storage savings compared to the original dense matrix. The 
degraded performance is due to that the fact that each block will require a 
large smooth component as well as detail components which are a result of 
the artificial boundaries created by the block method. It should be stressed 
that in spite of these artificial high components, fidelity is maintained as 
shown in figure 9. These observations are shown by comparing figure 10 
and 11 which is the impedance matrix of a 4A flat strip when the DAUB 10 
wavelet transform is applied on the entire matrix and in a block-by-block 
fashion with sixteen 32x32 blocks, respectively. 
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4 Integral Equations 

The matrix compression technique described in the previous section is a 
rather indirect use of wavelets in electromagnetics. In this section, we will 
look at a more direct application of wavelets for the solution of integral 
equations. This method will not prove to be very useful but will provide a 
reference for future research. 

The E-polarized scattering by a metallic strip may be solved by deter- 
mining the current via the integral equation 

£;<*, ») = j j c + (» - Vy\)d( (13) 

where the integration domain, C, suggests the parameterization 

x = x(t) ( 14 ) 

y = y(0 

Thus (13) may be written 

£(<) = 7 1 /*(‘>J , ’[*.|r(i,(’)|]<ft' (is) 

and 

r (M) = v^(0-*(O) 2 + (y(0-y(O) 2 (16) 

Let us now employ the method of weighted residuals using the discrete 
wavelet as an entire domain weight 

j C ' E\{t)xl> 2 ,(t - m2->)dt = jJ c J c Jz{t)Hi l \k 0 \r(tj‘)\\dt^ 2] (t-m2- 3 )dt{\l) 

If we now' expand the unknown current in terms of wavelet basis(8) and 
realize that all of the resulting integrals are DWTs, we arrive at the matrix 
equation 

ZJ* = £ (18) 

where now' the impedance matrix is given directly by a two-dimensional DWT 
of the Hankel function of (17) 

2 = VlW[*„|r(f,f')||] (19) 
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and the unknown vector is composed of modal coefficients. This solution 
is analogous to finding a Fourier series expansion solution if the complex 
exponential is used rather than the wavelet. We also note that the scattered 
field is computed efficiently as 


Kk 0 p 


N-l 




n=r 0 




( 20 ) 


where the excitation is computed via the FWT 

£x _ yy r e -'Meo5<J>x(()+s>n«!>y(< 


( 21 ) 


Figure 12 illustrates that the aforementioned formulation does recover the 
solution for a 1A metallic strip when 64 coefficients are determined. This 
is an unfortunately large number of terms required. The usefulness of this 
method does not improve with problem size as seen in figure 13 where even 
1024 coefficients are not sufficient to characterize a 150A strip. The failings 
of this particular formulation will be discussed in the next section. 


5 Conclusions and Future Work 

This paper has presented several different uses for wavelets or families of 
basis functions for computational electromagnetics. The relevant features 
and characteristics of wavelets were briefly discussed. Since the scope of 
this paper precludes a detailed review of wavelets, we mention that Strang[9] 
and Daubechies[10] give excellent reviews of wavelets. Mallat[l] not only 
introduces a particularly useful interpretation of wavelets and the wavelet 
transform, he developed from this viewpoint a quite efficient Fast Wavelet 
Transform which rivals the FFT in speed. Mallat also served to present this 
field in terms which are more accessible to engineers as compared the 
mathematicians “flavor 71 given by Daubechies. 

We developed a novel technique adapted from image processing for gener- 
ating a sparse matrix where previously a dense matrix was required. It should 
be stressed that the utility of this method solely is based of the scale charac- 
teristics of the impedance matrix rather than some restrictive mathematical 
operation such as convolution. The block-by-block compression method pro- 
posed resulted in significant storage savings. In the future, we are interested 
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in improving this method with respect to the degraded performance caused 
by block processing. It is reasonable to anticipate some improvement since 
only the most readily implemented techniques and analyzing wavelets were 
employed for this study and we understand the cause of the degradation. 
However it should be pointed out that if we can approach the efficiency of 
the entire matrix method, the block by block technique could be used with 
any existing code and a sparse matrix solver to dramatically extend the size of 
problems which could be handled with integral equation methods regardless 
of the geometry shape. 

The second method looked at was a direct solution of the integral equa- 
tion by wavelet expansion. This proved very disappointing. It is known 
that a Fourier series expansion of the current requires the retention of many 
terms. However, it was felt that the compact support property of wavelets 
would result in a dramatic reduction in the number of required coefficients. 
Although this initial trial did not succeed, we note that the wavelets cho- 
sen were available rather than the best possible. Further studv may reveal 
more appropriate analyzing wavelets for our purpose or perhaps higher order 
wavelets are required for sufficient approximation. 

Another avenue of research is employing wavelet testing with sub-domain 
basis functions. One may think of this operation as row-only compression 
rather than row and column compression which would result from the two- 
dimensional processing as described above. However, there would be no 
artificial components along the rows and by inspecting figure 11 we note 
that a significant amount of the spurious components are along rows of the 
matrix. This should result in a sparse matrix due to the local support of the 
basis and the compact support of the weights and a particularly appealing 
choice of basis are the ones used by Canningfll] since this basis will aid in 
compression along columns of the resulting matrix. To this date, we have not 
investigated this method. In addition, further study of Alpert et <z/[ 7] may 
yield further insight. It is noted that [7] investigated non-oscillating kernel 
functions while our kernel is certainly still a matter of research. 

Wavelet processing and direct solution offers a promising avenue for ex- 
tending integral equation techniques to much larger structures. It is noted 
that all methods described herein used FWT to evaluate all the integrals in 
an efficient manner. Although the examples presented involved the MoM, we 
are interested in using these techniques when terminating a Finite Element 
mesh with a boundary integral. An efficient solution in terms of unknowns 
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may be found for arbitrarily closures. This is in contrast the few available 
closures if one employs the CGFFT method to solve the boundary integral. 
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Table 1 


lambda Corner 

Reflector 

using 32 Unknowns 

and tol=10 A -2 


Wavelet 

Zeros 

Non-zeros 

Total 

Order 

DAUB4 

370 

654 

1 024 

20N 

DAUBG 

451 

573 

1 024 

20N 

DAUB8 

454 

570 

1 024 

20N 

DAUB10 

444 

580 

1 024 

20N 

DAUB12 

439 

585 

1024 

20N 

DAUB14 

436 

588 

1 024 

20 N 

DAUB16 

426 

598 

1024 

20N 

DAUB18 

420 

604 

1 024 

20N 

DAUB20 

387 

637 

1 024 

20N 


Table 2 


lambda Curved 

Strip using 

1 28 unknowns 

and tol=10' v -2 


Wavelet 

Zeros 

Non-zeros 

T otal 

Order 

DAUB4 

10311 

6073 

16384 

50N 

DAUB6 

12312 

4072 

1 6384 

30N 

DAUB8 

13214 

3170 

1 6384 

20N 

DAUB10 

1 3388 

2996 

16384 

20N 

DAUB12 

1 3597 

2787 

16384 

20N 

DAUB14 

1 3591 

2793 

16384 

20N 

DAUB16 

1 3643 

2741 

1 6384 

20N 

DAUB18 

1 3713 

2671 

16384 

20N 

DAUB20 

13656 

2728 

1 6384 

20N 



Table 3 

Resislive Junction using 1024 unknowns and tol=10 A -4 


Wavelet 

Zeros 

Non-zeros 

Total 

Order 

DAUB4 

573834 

474742 

1048576 

460N 

DAUB6 

777806 

270770 

1048576 

260N 

DAUB8 

871382 

177194 

1048576 

170N 

DAUB10 

917716 

130860 

1.048576 

130N 

DAUB12 

940642 

107934 

1048576 

1 10N 

DAUB14 

954802 

93774 

1 048576 

90N 

DAUB16 

966321 

82255 

1 048576 

80N 

DAUB18 

973085 

75491 

1048576 

70 N 

DAUB20 

978972 

69604 

1048576 

70N 


Resistive Junction using 
Wavelet Zeros 

Tabel 4 

block method w/ 
Non-zeros 

block size=64x64 
Total 

Order 

DAUB4 

268401 

7801 75 

1048576 

760N 

DAUB6 

432037 

616539 

1048576 

600N 

DAUB8 

569457 

479119 

1048576 

470N 

DAUB10 

685377 

363199 

1048576 

350N 

DAUB12 

744733 

303843 

1048576 

300N 

DAUB14 

786832 

261744 

1048576 

260N 

DAUB16 

796331 

252245 

1048576 

250N 

DAUB18 

797724 

250852 

1048576 

240N 

DAUB20 

799690 

248886 

1048576 

240N 
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Figure 3. a) h 0 ( x) - X| - i/2.i/n( x )» 

b) h x {x) - 2lH-(n)x lw/ 3.i/4.«/2^i/-!( x )- 

c) h x is decomposed into its "components", each component is a multiple of the characteristic 
function of an interval of length (The "components" of h j would have width 2"') 

d) Each "component" is replaced by a proportional version of centered around the same 
point as the component, and scaled down by a factor } (This scaling factor would be 2~ for h / * 

e) The functions in d) are added to constitute h z (or h t , if one starts from h } in c)) 


1/.x1a L-shaped Comer Reflector Impedance Matrix before Wavelet Transform 
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lhU L-shaped Comer Reflector Impedance Matrix after DAUB10 Wavelet Transform 
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